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ABSTRACT 
In this paper, the VSS T? control chart for monitoring the process mean is economically designed. The cost model proposed by 
Costa and Rahim (2001) to economically design the variable parameter X control chart is used here and it is minimized through a 
genetic algorithm (GA) approach. Then a comparison between the FRS and VSS schemes are made to clarify the economic 
advantages of the VSS scheme. 


Keywords: Hotelling’s T* control chart, Variable Sample size (VSS) scheme, Economic Design, Markov chain approach, Genetic 
Algorithm (GA). 


1. INTRODUCTION 

Nowadays there is growing interest in multivariate statistical quality control, i.e. the simultaneous control of several related quality 
characteristics of a process, because quality control problems in industry may involve more than a single quality characteristic, i.e. a 
vector of characteristics. This has formed the basis of extensive work performed in the field of multivariate quality control. Shewhart, 
who is famous for the development of the statistical control chart (Shewhart charts) first recognized the need to consider quality 
control problems as multivariate in character. A great deal of work on multivariate statistical control procedures was performed in 
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the 1930's and in the 1940's by Hotelling (1947). Jackson (1985) mentioned in his paper that the multivariate techniques should 
posse three important properties: (1) they produce a single answer to the question: is the process is in-control? , (2) has the specified 


type | error been mentioned?, and (3) these techniques must take into account the relationship between the variables. The 
Hotelling’s T* chart, satisfies the above three properties and also has the advantage of its simplicity. Consider p correlated 
characteristics are being measured simultaneously and are being controlled jointly. It is assumed that the joint probability 
distribution of the p quality characteristics is the p-variate normal distribution. The procedure requires computing the samples mean 


for each of the p quality characteristics from a sample of size 1. The subgroup statistics T* = n,(X—p,) > '(*- py) are plotted 


on a control chart in sequential order. Here, is the vector of in-control means for the p quality characteristics. For the sake of 


simplicity, in this paper, we assume here that the process mean vector p, and variance-covariance matrix >» are known. The 
statistic T? is distributed as a chi-square variable with p degrees of freedom. The upper limit on the control chart is UCL = x?(p)- 


When the process is in control, with a probability @ , the statistic T? exceeds the UCL, so that the overall error rate (type I) can be 
maintained exactly at the level @. If the process is out of control, the chart statistic is distributed as a non-central chi-square 


distribution with p degrees of freedom with non-centrality parameter 4 = n,d*, where d is the Mahalanobis distance that is used as 


a measure of process shift in multivariate statistical quality control. It is usually assumed that the variance-covariance structure of the 
quality characteristics being charted does not change and that the assignable cause is manifested by a shift or drift in at least one 
component of the mean vector of the process. 


The magnitude of this shift is often expressed by d* =(—p))' © ‘(uw—py), and p is the p-characteristics mean vector 7. The 


traditional practice in applying a control chart to monitor a process is to obtain samples of fixed size, 1). Therefore, traditional 


control charts use a Fixed Ratio Sampling (FRS) scheme. Variable Sample Size procedure (VSS) is a scheme that varies the sampling 
size as a function of prior sample results. The design of univariate Shewhart charts with adaptive sample sizes studied by Burr (1969), 
Daudin (1992), Prabhu et al. (1993), and Costa (1994). In these studies, two sample sizes are used and this procedure shows better 
power in detecting shifts in the mean. Zimmer et al. (1998) present a three-state adaptive sample size control chart and compared it 
with both the standard Shewhart control chart and the developed two-state adaptive sample size control chart. They conclude that 
three-state procedure is only slightly better than the two-state scheme, and so the two-state scheme is likely adequate in most 
applications 7. In the Multivariate SPC, The FRS T* control chart shows a good performance to detect large shifts in the process 
mean. However, in many practical situations it is necessary to be able to detect even moderate shifts in the process mean. In such 
cases, the statistical efficiency of the FRS T? chart (in terms of the speed with which process mean shifts are detected) is poor. Aparisi 
(1996) studied the VSS T? control chart. He considered an adaptive strategy for the subgroup size based on the data trends. He 
divided the area between the UCL and the origin, into two zones for the use of two different sample sizes. If the current sample 
value falls in a particular zone, then the corresponding sample size is to be applied for the successive sampling. He showed that the 
Hotelling’s T? control chart with VSS scheme, significantly improves the efficiency of the standard Hotelling’s T? control chart in 
detecting small changes in the process mean. For moderate shifts the obtained Average Run Length (ARL) values are about 1/3 of 
FRS scheme. However, Aparisi (1996) found that, when the shifts are large, ARL for the FRS chart is smaller or almost equal for the 
VSS chart 7. 

Aparisi (1996) assumed that the process goes out of control in the beginning of the monitoring the process. Faraz and 
Moghadam (2008) extended the work of Aparisi (1996) to the case in which the shift in the process mean does not occur at the 
beginning but at some random time in the future. Further, they assumed that the occurrence time of the shift is an exponentially 
distributed random variable. They showed that the proposed modification can improve the power of the chart to detect small to 
specifically moderate shifts in the process mean. Recently the economic design (ED) of control charts is of a great importance. Based 
on the ED procedure, the charts is designed in such a way that the overall costs associated with maintaining current control of a 
process is minimized. Weiler (1952) was one of the first researchers who attempted to design a control chart based on an economic 
viewpoint. He suggested that the optimal sample size should be the amount that minimizes the total number of inspected items 
required to detect a shift in the process mean. Four years later, Duncan (1956) proposed a very popular procedure for the ED of 
Shewhart X chart. Lorenzen and Vance (1986) developed an economic model for all kinds of control charts, not just the X chart. 
These economic models consider classical sampling scheme (FRS) and adopting these models to variable ratio sampling (VRS) 
procedures seems somewhat difficult. Costa and Rahim (2001) developed an economic model based on the Markov chain approach 
which is suitable to study the economic design of the VRS schemes. 

Montgomery (1980) made a very good literature review on ED of control charts and listed fifty one references on the topic. 
Woodall (1986) criticized the EDs by noting that these models result in designs with high Type | error. Saniga (1989) amended this 
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drawback by adding some statistical constraints on the Type | and/or Type Il errors. He called this approach the Economic Statistical 
Design (ESD). He observed that the ESD statistical properties are as good as statistically designed control charts. We hereby propose 
the economic design of the VSS T? control schemes which has not been found in the literature. The paper is organized as follows: In 
section 2 the VSS T* control scheme proposed by Faraz and Moghadam (2008) is briefly reviewed and upon the Markov chain 
approach, the cost model proposed by Costa and Rahim (2001) is constructed. Section 3 makes a comparison between the FRS and 
VSS schemes. Finally, concluding remarks make up the last section. 


2. THE COST MODEL 
Consider p correlated quality characteristics are to be monitored simultaneously using the T* scheme. First, it is assumed that the 
joint probability distribution of the quality characteristics is a p-variate normal distribution with the known in-control mean vector 


Hi) and known variance-covariance matrix >! . Then the subgroups’ statistic T, =n(X, —ji,)' © *(%, —ji,) (each of size n) is 
plotted in sequential order and the chart signals as soon as T? >k = y7(p). To gain from the Markov chain properties, following 


transient states are defined: 


State 1: 0 < T* < w and the process is in control; 
State 2: w<T? <kand the process is in control; 
State 3: T? >k and the process is in control; 

State 4: 0 < T* < w and the process is out of control; 
State 5: w<T? <k and the process is out of control; 


State 6: T? >k and the process is out of control. 


The state 3 stands to count the false alarms and the state 5 represents the absorbing state. Hereafter, the transition probability 
matrix is given by 

Pu Po Ps Pi Piss Pre 

Pox P22 P23 Po Pr Pr 

P P31 Px. P33 P34 P35 P36 

0 0 O Pus Pas Pas 


0 0 O Psy DPss Ps 
0) 0) 0 0 0 1 


where P;; denotes the transition probability that i is the prior state and j is the current state. In what follows, F(x, p,7) will 
denote the cumulative probability distribution function of a non-central Chi-square distribution with p degrees of freedom and non- 


centrality parameter 77 = n,d*, where d is the amount of the process mean shift presented by Mahalanobis distance measure. 


The P;;'s are given as follows: 


Ah 


1 = Px = Ps = Fw, p,0) xe 
2 = Px = Px = Fk, p,0)xe™ — p,, 
<= P4=04Se" —Pe= Pa 

4 = F(w, p,n,d’)x(l-e”) 

5 = F(k, p,n,d*)x(l-e™) - py, 


Ah 


gale =p Pe 


v DD DDS 
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Pu = Py = Fw, p.n,d?)xd—e*) 

Pos = Pas = F(k, p,n,d*)x(1-e") — py, 
Pog = Psp = 1-6 — Dos — Pog 

Pa = F(w, p,n,d’) 

Pas = F (Kk, p,m’) — Dag 

Pas =1— Pas — Pas 

Psa = F(w, p,n,d’) 

Pss = F (k, pnd”) — Psg 

Ps5 =1— Pss — Psa 


Now, upon the cost model specified by Costa and Rahim (2001), the expected loss per hour is determined by following function. 
In what follows, V, is the hourly profit earned when the process is in control, V, is the hourly profit earned when the process is out- 


of-control, Co is the expected cost for investigating false alarms, C7 is the expected cost to find assignable cause and repair the 
process, s stands for the cost of each inspected item, 7, determines the time needed to investigate a false alarm, T,is the time 


needed to search for and repair the assignable, b’=(0, 1, 0, 0, 0) is a vector of initial probabilities, I is the identity matrix of order 
5, Qis the 5x5 matrix obtained from P on deleting the elements corresponding to the absorbing state, h’=(h, h, h, h, h) is the 
vector of sampling intervals, i, =(0, 0, 1, 0, 0)’ and n=(m, nm, 1, m, 13)’ 


Y—V 


E(A) =V, 


+b'(I-Q)'(V,xh-C, x1, —sxii)—-C, (1) 


b'(I-Q) '(h+7, x1,) +7, 


In order to economically design the VSS, we try to find five chart parameters (control limit k, warning line W, sample sizes 1, 
and ny and sampling frequency) that minimize E(L). The cost function (1) has a discontinuous and non-convex solution space and 


hence is solved by using genetic algorithms (GAs). The crossover rate 0.1, the population size 100 and the mutation rate 0.5 are used 
here and the number of generation is set to at least 10,000 times. 


3. COMPARISON RESULTS 

Table 1 gives the thirteen process and cost parameters which is taken from Costa and Rahim (2001). These values provide a general 
variation in parameter values and shall be used here for our comparison purpose. Table 2 gives the optimal design parameters and 
resulting loss for the economic VSS scheme with a comparison to the economic FRS scheme for the cases p=2. The results are 
opposed to what Aparisi (1996) and Faraz and Moghadam (2008) achieved. They concluded that the VSS scheme detects the 
process mean shifts more quickly than the FRS scheme and hence reduce the resulting cost due to nonconforming items. Their 
result is somehow true but we show that when we take a look at the VSS scheme by considering complete economic information, it 
is not economical than the FRS scheme. In fact, the results indicate that the VSS scheme can be designed in such a way that is able 
to detect the process mean shifts more quickly than the FRS scheme, but clearly this method is not economical at all. Hence, the 
other VRS schemes like the VSI and the VSSVSI can be applied to achieve this goal. 


Table 1 The 13 process and cost parameters 


No. Ss Ci Cc Yo V, i Ee. A d 
5 | 500 | 500 |} 500 ; 50 5 1 0.01 1 

2 10 | 500 | 500 | 500 |} 50 5 1 0.01 1 
5 | 250 500 | 500 | 50 5 1 0.01 1 
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4 5 | 500 | 50 | 500 | 50 5 1 0.01 1 
5 5 | 500 | 500 | 250 | 50 5 1 0.01 1 
6 5 | 500 | 500 | 500 | 100 5 1 0.01 1 
7 5 | 500 | 500 | 500 0 5 1 0.01 1 
8 5 ; 500 | 500 | 500 | 50 | 2.5 1 0.01 1 
9 5 | 500 | 500 | 500 | 50 5 10 | 0.01 1 
10 5 | 500 | 500 | 500 | 50 5 1 0.05 1 
11 5 | 500 | 500 | 500 | 50 5 1 0.01 | 1.5 
12 5 | 500 | 500 | 500 | 50 5 1 0.01 | 0.5 
13 5 | 500 | 500 | 500 | 50 5 1 0.01 | 2.0 


Table 2 The optimal parameters of the economic design of the FRS and VSS schemes for the case p =2 


No FRS scheme VSS scheme 
k h n | ANF | AATS E(L) k w h n; | n2| ANF | AATS | E(L) 
1 | 10.49 | 6.27 | 18} 0.08 | 4.08 | 43.37 | 10.46 | 4.96 | 5.86 | 16 |] 20 | 0.08 | 4.05 | 43.23 
2 9 8.37 | 16] 0.13 | 545 | 5462 | 9.08 | 439 | 812 | 15] 18] 0.12 | 5.42 | 54.49 
3 | 10.31} 622 | 18] 0.09 | 4.00 | 43.17 | 10.44 | 5.22 | 6.04 | 17 | 21 | 0.09 | 3.98 | 43.05 
4 | 10.51 | 624 | 18} 0.08 | 4.07 | 39.10 | 10.64 | 5.12 | 6.06 | 17 | 21 | 0.08 | 4.05 | 38.97 
5 9.32 | 9.12 | 16] 0.10 | 6.05 | 2854 | 942 | 42 | 887 | 15 | 18} 0.1 6.01 | 28.48 
6 10.5 | 668 | 18 | 0.08 | 4.35 | 41.38 | 10.64 | 5.14 | 649 | 17] 21] 0.08 | 4.32 | 41.27 
7 10.5 | 5.94 | 15 | 0.09 | 3.87 | 45.25 | 10.64 | 5.14 | 5.76 | 14] 18] 0.09 | 3.84 | 45.11 
8 | 943 | 5.97 | 16] 0.14 | 4.00 | 42.22 | 949 | 4.72 | 5.78 | 15 | 18 | 0.14 | 3.98 | 42.1 
9 | 10.36 | 6.53 | 18} 0.08 | 4.22 | 79.24 | 10.49] 5.2 | 633 | 17] 21] 007 | 42 79.13 
10 | 10.03 3 17 | 0.04 ) 2.01 | 114.23 | 10.2 | 4.68 | 3.03 | 17 | 19 | 0.04 2 114.1 
11 | 12.22 | 444 | 9 | 0.05 | 2.89 | 33.52 | 12.49 | 6.02 | 446 | 9 | 10 | 0.04 2.9 33.38 
12 | 7.46 | 11.06 | 52) 0.21 | 7.59 | 68.95 | 7.46 | 0.00 | 11.06 | 25 | 52 | 0.21 | 7.59 | 68.95 
13 | 13.42 | 3.46 | 5 | 0.03 | 2.07 | 2830 | 13.62 | 664] 3.47 | 5 | 6 | 0.02 | 2.06 | 28.24 


4. CONCLUDING REMARKS 
The Hotelling's T? control chart is a widely applied multivariate control chart for detecting moderate to large shifts in process mean. 


In the present paper, the economic design of the VSS T* control chart is developed based on the cost model considered by Costa 


and Rahim (2001). The expected total cost per hour is minimized using GA and finally a comparison is made between the FRS and 


VSS schemes from an economic viewpoint. The results indicate that the VSS scheme is not more economical than the FRS scheme 


and hence the VSS scheme is not preferred to the FRS scheme for monitoring purpose. 
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